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Abstract 

We present realistic models for massive Population III 
stars and stellar populations based on non-LTE model at- 
mospheres, recent stellar evolution tracks and up-to-date 
evolutionary synthesis models, with the aim to study their 
spectral properties, including their dependence on age, 
star formation history, and IMF. 

A comparison of plane parallel non-LTE model at- 
mospheres and comoving frame calculations shows that 
even in the presence of some putative weak mass loss, the 
ionising spectra of metal-free populations differ little or 
negligibly from those obtained using plane parallel non- 
LTE models. As already discussed by Tumlinson & Shull 
(2000), the main salient property of Pop III stars is their 
increased ionising flux, especially in the He"*" continuum 
(> 54 eV). 

The main result obtained for individual Pop III stars is 
the following: Due to their redward evolution off the zero 
age main sequence (ZAMS) the spectral hardness mea- 
sured by the He+/H ionising flux is decreased by a factor 
^ 2 when averaged over their lifetime. If such stars would 
suffer strong mass loss, their spectral appearance could, 
however, remain similar to that of their ZAMS position. 

The main results regarding integrated stellar popula- 
tions are: 

— For young bursts and the case of a constant SFR, neb- 
ular continuous emission neglected in previous stud- 
ies — dominates the spectrum redward of Lyman-a if 
the escape fraction of ionising photons out of the con- 
sidered region is small or negligible. In consequence 
predicted emission line equivalent widths are consid- 
erably smaller than found in earlier studies, whereas 
the detection of the continuum is eased. Nebular line 
and continuous emission strongly affect the broad band 
photometric properties of Pop III objects. 

— Due to the redward stellar evolution and short lifetimes 
of the most massive stars, the hardness of the ionising 
spectrum decreases rapidly, leading to the disappear- 
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ance of the characteristic He ii recombination lines af- 
ter ~ 3 Myr in instantaneous bursts. 

— He II A1640, Ha (and other) line luminosities usable 
as indicators of the star formation rate are given for 
the case of a constant SFR. For obvious reasons such 
indicators depend strongly on the IMF. 

— Due to an increased photon production and reduced 
metal yields, the relative efficiency of ionising photon 
energy to heavy element rest mass production, ry, of 
metal-poor and metal-free populations is increased by 
factors of ^ 4 to 18 with respect to solar metallicity 
and for "standard" IMFs. 

— The lowest values of 77 ^1.6 2.2% are obtained 
for IMFs exclusively populated with high mass stars 
(Miow 50 Mq). If correct, the yields dominated by 
pair creation SNae then predict large overabundances 
of 0/C and Si/C compared to solar abundance ratios. 

Detailed results are given in tabular form and as fit for- 
mulae for easy implementation in other calculations. The 
predicted spectra will be used to study the detectability 
of Pop III galaxies and to derive optimal search strategies 
for such objects. 

1. Introduction 

Important advances have been made in recent years on 
the modeling of the first stars and galaxies forming out 
of pristine matter — so called Population III (Pop III) 
objects — in the early Universe (see e.g. the proceedings 
of Weiss et al. 2000 and Umemura & Susa 2001). Among 
the questions addressed (see also the reviews of Loeb & 
Barkana 2001 and Barkana & Loeb 2001) are the star 
formation process and the initial mass function (IMF) of 
the first stars (Tegmark et al. 1997, Abel et al. 1998, 2000, 
Bromm et al. 1999, Nakamura & Umemura 1999, 2001, 
Omukai & Palla 2001), their effect of energy injection and 
supernovae (SN) (MacLow & Ferrara 1999, Ciardi et al. 
2000), their role on cosmic reionisation (Haiman & Loeb 
1997, Gnedin & Ostriker 1997, Tumlinson & Shull 2000, 
Cojazzi et al. 2000, Ciardi et al. 2001), metal enrichment of 
the IGM and other signatures of early chemical evolution 
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(Ferrara et al. 2000, Abia et al. 2001), and dust formation 
(Todini & Ferrara 2000). 

Extensive sets of recent stellar evolution and nucle- 
osynthesis calculations for zero metallicity stars covering 
a wide range of stellar masses have also become available 
recently (Marigo et al. 2000, Feijoo 1999, Desjacques 2000, 
Woosley & Weaver 1995, Heger et al. 2000, 2001, Umeda 
et al. 2000). 

Tumlinson & ShuU (2000, hereafter TSOO) have re- 
cently pointed out the exceptionally strong He"*" ionising 
flux of massive (M ^ 40 Mq) Pop III stars, which must 
be a natural consequence of their compactness, i.e. high 
effective temperatures, and non-LTE effects in their atmo- 
spheres increasing the flux in the ionising continua. As a 
consequence strong He ii recombination lines such as He ii 
A1640 or He ii A4686 are expected; together with AGN a 
rather unique feature of metal-free stellar populations, as 
discussed by TSOO and Tumlinson et al. (2001, hereafter 
TGSOl). Instead of assuming a "standard" Population I 
like Salpeter IMF and "normal" stellar masses up to 100 
Mq as TSOO, Bromm et al. (2001) have considered stars 
with masses larger than 300 Mq , which may form accord- 
ing to some recent hydrodynamical models (Abel et al. 
1998, Bromm et al. 1999, Nakamura & Umemura 2001). 
An even stronger ionising flux and stronger H and He ii 
emission lines was found. 

Some strong simplifying assumptions are, however, 
made in the calculations of TSOO, TGSOl, and Bromm 
et al. (2001). 

1) All stars are assumed to be on the zero age main se- 
quence, i.e. stellar evolution proceeding generally to 
cooler temperatures is neglected 0. 

2) Bromm et al. (2001) take the hardness of the stellar 
spectrum of the hottest (1000 Mq) star as represen- 
tative for all stars with masses down to 300 Mq. This 
leads in particular to an overestimate of the He ii re- 
combination line luminosities. 

3) None of the studies include nebular continuous emis- 
sion, which cannot be neglected for metal-poor objects 
with such strong ionising fluxes. This process increases 
significantly the total continuum flux at wavelengths 
redward of Lyman-a and leads in turn to reduced emis- 
sion line equivalent widths. 

4) A single, fixed, IMF is considered only by TSOO and 
TGSOl. In view of the uncertainties on this quantity 
it appears useful to explore a wider range of IMFs. 

The models presented here relax all the above as- 
sumptions and allow the exploration of a wide parameter 
space in terms of stellar tracks (including also alternate 
tracks with strong mass loss), IMF, and star formation 
history (instantaneous bursts or constant star formation 
rate (SFR)). For this aim we provide in particular the 
main results in tabulated form and as fit formulae. 

^ A simple estimate of this effect is given by TGSOl. 



First we examine the influence of several poorly con- 
strained physical processes (mass loss, non coherent elec- 
tron scattering etc.) on the predicted stellar fluxes of Pop 
HI stars using non-LTE atmosphere models. Combined 
with metal-free stellar tracks these predictions are then 
introduced to our evolutionary synthesis code (Schaerer 
& Vacca 1998) updated to calculate the nebular proper- 
ties (H and He emission lines and nebular continuum) for 
metal-free gas. In this manner we examine the proper- 
ties of both individual Pop III stars and integrated stellar 
populations. 

Calculations of the relative photon to metal produc- 
tion of stellar populations (e.g. Madau & Shull 1996) have 
proven useful for a variety of studies ranging from esti- 
mates of the stellar contribution to the UV background 
radiation (Madau & Shull 1996), over studies of metals 
in the IGM, to SN rates in the early Universe and their 
detectability (Miralda-Escude & Rees 1997). To improve 
on such calculations based on solar metallicity models we 
also calculate the heavy element production from metal- 
poor and metal-free populations including contributions 
from type II SN and possible pair creation SN from very 
massive stars (cf. Ober et al. 1983, Heger et al. 2000). 

It is the hope that our models should in particular help 
to examine more precisely the properties of the first lumi- 
nous objects in the Universe in preparation of their first 
direct observation (cf. TGSOl, Oh et al. 2001a, Schaerer 
& Pello 2001, PcUo & Schaerer 2001). 

The paper is structured as follows: The main model 
ingredients are summarised in Sect. ^ In Sect. ^ we dis- 
cuss the properties of individual Pop III stars. The proper- 
ties of integrated stellar populations, their spectra, photon 
fluxes, metal production etc. are presented in Sect. Our 
main results are summarised in Sect. ^. 

2. Model ingredients 

2.1. Atmosphere models 

As well known (e.g. Mihalas 1978) strong departures from 
local thermodynamic equilibrium (LTE) occur in the at- 
mospheres of hot stars. In addition, for the most massive 
stars it is of interest to examine the possible effect of mass 
loss on their emergent spectra. Albeit presumably small, 
mass loss — due to radiation pressure or pulsational insta- 
bilities (Kudritzki 2000, Baraffe et al. 2000) — potentially 
affects the spectral energy distribution (SED), especially 
in the ionising continua, as known for Pop I stars (Gabler 
et al. 1989, Schaerer & de Koter 1997). For stars with 
negligible ionising fluxes, the remainder of the SED of in- 
terest here (rest-frame UV and optical) is well described 
by plane parallel LTE atmospheres. 

To account for these cases we use the following model 
atmospheres: 

— For stars with T^ff > 20000 K we constructed an ex- 
tensive grid of pure H and He plane parallel model at- 
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mospheres using the non-LTE code TL US TY code of 
Hubeny & Lanz (1995). As far as it is possible to con- 
struct hydrostatic structures, we calculated a model 
grid with various \ogg values (typically Alogg = 0.5) 
to cover the domain corresponding to the evolutionary 
tracks (see below). The primordial He abundance of 
Izotov & Thuan (1998, mass fraction Y = 0.244) was 
adopted. The influence of the exact H/He abundance 
on the spectra is small. 

— To explore the possible importance of mass loss on the 
ionising spectra of massive Pop III stars, we use the 
CMFGEN code of Hillier & Miller (1998), which calcu- 
lates the radiation transfer in the co-moving frame and 
solves for non-LTE equilibrium in radiative equilib- 
rium. The input atmospheric structure including pho- 
tosphere and wind is calculated with the ISA- WIND 
code of de Koter et al. (1996). Models were calculated 
for stellar parameters corresponding to the zero age 
main sequence (ZAMS) of 25, 60, and 500 Mq stars 
(parameters given in Sect. |3.2| ) and an evolved 60 M© 
model at 52 kK, and wind parameters foo ~ 2000 km 
s-i and M = 10-9 to 10"^ Mq yr"! (10"^ to 10"^ 
for the 500 Mq model). The same H/He abundances 
as above are adopted. 

For stars with strong mass loss (as also explored below) 
we have used the pure Helium Wolf-Rayet (WR) atmo- 
sphere models of Schmutz et al. (1992). The two pa- 
rameters, core temperature and wind density, required 
to couple these models to the stellar evolution models 
are calculated as in Schaerer & Vacca (1998). 

— The plane parallel LTE models of Kurucz (1991) with 
a very metal-poor composition ([Fe/H] ~ ~5.) are 
adopted for stars with T^s < 20000 K. The models 
with such a low metal content are equivalent to pure 
H and He LTE models (Cojazzi et al. 2000, Chavez et 
al. 2001). 
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Fig. 1. HR-diagram for metal free {Z = 0, solid and long- 
dashed lines) and low metallicity (Z — 1/50^©, dotted) 
stars. Isochrones of 2 and 4 Myr for Z — tracks with- 
out mass loss are also plotted (long-dashed). The short 
dashed high mass tracks evolving blueward of the ZAMS 
are computed assuming strong mass loss. The position of 
the ZAMS at solar metalhcity (Schaller et al. 1992) is 
shown by the solid line. Pop III tracks and isochrones from 
Marigo et al. (2001, no mass loss) and Klapp (1983), and 
El Eid et al. (1983, strong mass loss). Z = 0.0004 tracks 
for 0.8 to 150 Mq from Lejeune & Schaerer (2001). Note 
the important shift of the ZAMS to high Toff from low 
metallicity to Z = 0, as well as the rapid redward evolu- 
tion of the massive stars 



2.2. Stellar tracks 

To explore two extreme evolutionary scenarios for massive 
stars we have constructed two different sets of zero metal- 
licity stellar evolution tracks: 1) no mass loss, 2) strong 
mass loss (see Fig. |l|). 

For tracks with no or negligible mass loss (1) we 
use recent tracks from 1 to 500 Mq calculated with the 
Geneva stellar evolution code (Feijoo 1999, Desjacques 
2000). These tracks, including only the core H-burning 
phase and assuming small mass loss, have been compared 
with the tracks of Marigo et al. (2001) up to 100 Mq and 
additional models up to 500 Mq by Marigo (2000, pri- 
vate communication). Good agreement is found regard- 
ing the zero age main sequences (ZAMS), H-burning life- 
times, and the overall appearance of the tracks. The dif- 
ferences due to mass loss in the Geneva models are mi- 
nor for the purpose of the present work. Furthermore, as 



the He-burning lifetime is ^ 10 % of the main sequence 
phase and is spent at cooler temperatures (cf. Marigo et al. 
2001), neglecting this phase has no consequences for our 
predictions. This has also been verified by comparisons of 
integrated stellar populations adopting alternatively the 
isochrones provided by Marigo et al. (2001). 

The "strong mass loss" set (2) consists of a combina- 
tion of the following stellar tracks for the mass range from 
80 to 1000 Mq. Tracks of Klapp (1983) for initial masses 
of 1000 and 500 Mq computed with values of iV = 50 
and 100 respectively for the mass loss parameter^. For 
300, 220, 200, 150, 100, and 80 M© we adopt the tracks 
of El Eid et al. (1983). The remaining models (1 to 60 
Mq) are from the "no mass loss" set. The main differ- 
ence between the "strong" and "no mass loss" sets is the 
rapid blueward evolution of the stars in the former case, 

^ Tracks with these A'' values yield a qualitatively similar 
behaviour as the El Eid et al. (1983) models. 
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due to strong increase of the He abundance on the surface 
of these stars leading to a hot WR-like phase (see Fig. 
|l|). While the use of updated input physics (e.g. nuclear 
reaction rates) could lead to somewhat different results if 
recomputed with more modern codes, the predicted tracks 
depend essentially only on the adopted mass loss (see also 
Chiosi & Maeder 1986), and remain thus completely valid. 

In Figure |l| the two set of tracks are compared to the 
low metallicity (Z = 1/50 Z0, dotted lines) models of 
Lejeune & Schaerer (2000) and the position of a solar 
metallicity ZAMS (Schaller et al. 1992). This illustrates 
the well known fact that - due to the lack of CNO ele- 
ments - the ZAMS of massive Pop III stars is much hotter 
that their solar or low metallicity counterparts (cf. Ezer 
& Cameron 1971, El Bid et al. 1983, Tumhson & ShuU 
2000). In particular this implies that at Z — stars with 
M <; 5 Mq have unusually high temperatures and in turn 
non- negligible ionising fluxes (cf. Sect. ||) corresponding 
to normal 0-type stars (Toff ^ 30 kK). 



2.3. Evolutionary synthesis models 

To calculate the properties of integrated zero metallic- 
ity stellar populations we have included the stellar atmo- 
sphere models and evolutionary tracks in the evolutionary 
synthesis code of Schaerer & Vacca (1998). The follow- 
ing changes have been made to adapt the calculations to 
metal- free populations. 



2.3.1. Nebular emission 

The potentially strongest H and He recombination lines in 
the rest-frame UV from Lyman-a longward to Ha in the 
optical spectrum are calculated, namely Lyman-a, He 11 
A1640, He 11 A3203, He i A4026, He i A4471, He 11 A4686, 
H/3, He I A5016, He i A5876, and Ha|. We assume Case B 
recombination for an electron temperature of Tg — 30000 
K appropriate of metal-free gas and a low electron density 
(rie — 100 cm^"^). The emissivities and recombination co- 
efhcients of Storey & Hummer (1995) are used in general. 
For He i A4471 we take the data from Osterbrock (1989) 
for Tg — 20000 K. Lyman-a emission is computed assum- 
ing a fraction of 0.68 of photons converted to Lyman-a 
(Spitzer 1978). For these assumptions the line luminosi- 
ties (in units of erg s~^) are thus expressed as 
Li [ergs-i] =Q (l-/csc) Q^ [s-'l (1) 
where c; are the coefficients given in Table 0, /esc is the 
photon escape fraction out of the idealised region consid- 
ered here (we assume fcsc = in our calculations, cf. Sect. 
^), and Qi is the ionising photon flux (in units of s^^) cor- 
responding to the appropriate recombination line. The line 
emissivities depend only relatively weakly on the assumed 
conditions. Adopting, e.g. a temperature of Tg — 10000 



Table 1. Line emission coefficients c; for Case B, Tg 
30000 K, ng = 100 cm'^ 



Line 


ci [org] 




appropriate ionising flux Qi 


Lyman- Of 


1.04 X 10" 


11 


0(H) 


He II A1640 


5.67 X 10" 


12 


Q(He+) 


He II A3203 


3.32 X 10" 


13 


Q(He+) 


He I A4026 


9.25 X 10~ 


13 


Q(He») 


He I A4471 


1.90 X 10~ 


13 


Q(He») 


He II A4686 


7.13 X 10" 


13 


Q(He+) 


H/3 


4.47 X 10" 


13 


Q(H) 


He I A5016 


1.31 X 10" 


13 


g(He°) 


He I A5876 


4.90 X 10" 


13 


Q(He») 


Ha 


1.21 X 10" 


12 


Q(H) 



K leads to an increase of the He 11 line luminosities by ^ 
10 - 30 %, whereas changes of ;^ 10 % arc found for the 
H lines (ci64o = 6.40 x 10-12, CHa = 1.37 x lO^^^). This 
lower value of Tg is adopted for calculations at non-zero 
metallicities as in Schaerer & Vacca (1998) 

As we will be studying objects with strong ionising 
fluxes, nebular continuous emission must also be included. 
The monochromatic luminosity of the gas is given by 

(2) 



1, = ^2:^2^i(i-/_)Q(H) 

A^ as 

(e.g. Osterbrock 1989), where a^ is the case B recombi- 
nation coefficient for hydrogen. The continuous emission 
coefficient 7totah including free-free and free-bound emis- 
sion by H, neutral He and singly ionised He, as well as the 
two-photon continuum of hydrogen is given by 

n(Hc+) n(He++) 

7total = 7HI + 72q + 7HcI + 7HcII ' (3) 

The emission coefficients 7^ (in units of erg cm"^ s^^ 
Hz^^) are taken from the tables of AUer (1984) and Fer- 
land (1980) for wavelength below/above 1 /im respectively 
for an electron temperature of 20000 K 0. The primor- 
dial H/He abundance ratio of Izotov & Thuan (1998) is 
adopted. As a non-negligible He 11 ionising flux is obtained 
in the calculations under certain conditions (cf . below) , we 
also allow for possible continuous emission from the region 
with doubly ionised He. The relative contribution of the 
He++ and the He+ regions is to first order easily calculated 
from the relative size of the respective Stroemgren spheres 
and is given by (n(He++)/n(H+))/(n(He+)/n(H+)) ~ 

2. X (5(Hc+)/Q(H). As expected, the resulting contribu- 
tion from the He+ region is small, even in the case of the 
hardest spectra considered here. However, as shown later, 
the total nebular continuous emission is not negligible for 
Population HI objects even in the rest-UV range. 



^ More precisely: all He 11 lines with relative intensities 
Ix/I-ai3 ^ 0.4, He I with Ix/Uaii ^ 0.4, plus H/3 and Ha. 



^ For non-zero metallicities we adopt again Te — 10000 K as 
in Schaerer & Vacca (1998). 
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2.3.2. Initial mass function 

The question of the masses and mass distribution of metal- 
free stars has been studied since the 1960's (see the volume 
of Weiss et al. 2000 for recent results). 

While the outcome from various hydrodynamical mod- 
els and other studies still differ, there seems to be an over- 
all consensus that stars with unusually large masses (up 
to ~ 10^ Mq) may form, possibly even preferentially (e.g. 
Abel et al. 1998, Bromm et al. 1999, 2001, Nakamura & 
Umemura 2001). Uehara et al. (1996) and Nakamura & 
Umemura (1999, 2001) also find that the formation of 
stars with masses down to ~ 1 M0 is not excluded (cf. 
also Abel et al. 2001). The differences may be of various 
origins (adopted numerical scheme and resolution, treat- 
ment of radiation transfer and optically thick regions etc. ) . 

In view of our ignorance on this issue we adopt a va- 
riety of different upper and lower mass limits of the IMF 
assumed to be a powerlaw, with the aim of assessing their 
impact on the properties of integrated stellar populations. 
The main cases modeled here are summarised below in 
Table |[ Except if mentioned otherwise, the IMF slope is 
taken as the Salpeter value {a = 2.35) between the lower 
and upper mass cut-off values Mjow and Mup respectively. 
The consideration of stars even more massive than 1000 
Mq (or no mass loss models with > 500 M©) is limited 
by the availability of evolutionary tracks for such objects. 
The tabular data given below allows the calculation of in- 
tegrated populations with arbitrary IMFs for ZAMS pop- 
ulations and the case of a constant SFR . 

3. Properties of individual stars 

We now discuss the properties of the ionising spectra and 
their dependence on assumptions of model atmospheres. 
The ZAMS properties of individual stars as well as their 
average properties taken over their lifetime (i.e. including 
the effect of stellar evolution) are presented next. 

3.1. Ionising and H2 dissociating spectra 

The ionising photon fluxes qi and the hardness of the spec- 
trum (described by q{lle°)/q{H) and q{Re~^)/q{H)) pre- 
dicted by various atmosphere models as a function of Toff 
are shown in Fig. |^. The ionising photon flux per surface 
area qi, given in units of photons s^^ cni^^, is related to 
the stellar ionising photon flux Qi (in photons s^^) used 
later by 

poo p 

Q, = AirRlq, = AttRI / -^d,y, (4) 

where Hi, is the stellar radius, F^, the spectral flux distri- 
bution (in erg cm~^ s^^ Hz^^), and Vi the frequency of the 
appropriate ionisation potential. The photon flux (5(LW) 
in the Lyman- Werner band (11.2 to 13.6 eV) capable to 
dissociate H2 is also calculated. 

The solid lines in Fig. || give the predictions obtained 
from the plane parallel non-LTE TLUSTY models. As al- 




Teft [kK] 

Fig. 2. Ionising photon flux (top panel) and hardness of 
the ionising spectrum (middle and bottom panel) as a 
function of effective temperature for various atmosphere 
models. Solid lines connect plane parallel TLUSTY Pop 
III composition models of various Toff and log g. Triangles 
show calculations from the spherically expanding CMF- 
GEN models for Pop III stars. For comparison the pure 
H and He WR models of Schmutz et al. (1992) are also 
shown (open squares and circles). Discussion in text 



ready discussed by Tumlinson & ShuU (2000, TSOO), this 
shows in particular the high fraction of photons emitted 
with energies capable to ionise He i and He 11 for stars 
with Teff ^ 40 kK and ^ 80 kK respectively. 

As the stellar tracks with strong mass loss evolve from 
the ZAMS (at - 100 kK for M ^ lOOM©) blueward, we 
have also extended the temperature range to the hottest 
temperatures for which atmosphere models with winds are 
available (Schmutz et al. 1992, shown as open symbols in 
Fig. ph. As best shown in the bottom panel, the ratios 
(7(He")/g(H) and g(He+)/g(H) continue to increase with 
Tcff, as the maximum of the flux distribution has not yet 
shifted above 54 eV. However, it is important to remember 
that in the case of atmospheres with winds, the ionising 
flux also depends on the wind density as illustrated by the 
dispersion of the Schmutz et al. models for given T^s. For 
example, for sufflciently dense winds the He 11 ionising flux 
can be completed suppressed (cf. Schmutz et al. 1992). 
However, for the wind densities obtained in the present 
models, this situation is not encountered (cf. below). 
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Table 2. Summary of main model input parameters. All models assume a Salpeter slope for the IMF 
Model ID tracks Miow Mup symbol in Figs. |, 0, | 



A No mass loss 


1 


100 




dotted 


B No mass loss 


1 


500 




solid 


C No mass loss 


50 


500 


lonj 


5 dashed 


D Strong mass loss 


1 


500 


dash 


- dotted 


E Strong mass loss 


50 


1000 


long dash 


- dotted 


Solar metallicity tracks: 










ZS High mass loss 


1 


100 






Z = 1/50 Zq tracks: 










ZL Mass loss 


1 


150 







Table 3. Properties of Pop III ZAMS stars 



Mini 


log L/Lq 


logTeff 


Q(H) 


Q{Be°) 


Q(He+) 


Q(LW) 


Q(He'')/Q(H) 


Q(He+)/Q(H) 


1000. 


7.444 


5.026 


1.607E+51 


1.137E+51 


2.829E+50 


1.727E+51 


0.708E+00 


0.176E+00 


500. 


7.106 


5.029 


7.380E+50 


5.223E+50 


1.299E+50 


7.933E+50 


0.708E+00 


0.176E+00 


400. 


6.984 


5.028 


5.573E+50 


3.944E+50 


9.808E+49 


5.990E+50 


0.708E+00 


0.176E+00 


300. 


6.819 


5.007 


4.029E+50 


2.717E+50 


5.740E+49 


4.373E+50 


0.674E+00 


0.142E+00 


200. 


6.574 


4.999 


2.292E+50 


1.546E+50 


3.265E+49 


2.487E+50 


0.674E+00 


0.142E+00 


120. 


6.243 


4.981 


1.069E+50 


7.213E+49 


1.524E+49 


1.161E+50 


0.674E+00 


0.142E+00 


80. 


5.947 


4.970 


5.938E+49 


3.737E+49 


3.826E+48 


6.565E+49 


0.629E+00 


0.644E-01 


60. 


5.715 


4.943 


3.481E+49 


2.190E+49 


2.243E+48 


3.848E+49 


0.629E+00 


0.644E-01 


40. 


5.420 


4.900 


1.873E+49 


1.093E+49 


1.442E+47 


2.123E+49 


0.584E+00 


0.770E-02 


25. 


4.890 


4.850 


5.446E+48 


2.966E+48 


5.063E+44 


6.419E+48 


0.545E+00 


0.930E-04 


15. 


4.324 


4.759 


1.398E+48 


6.878E+47 


2.037E+43 


1.760E+48 


0.492E+00 


0.146E-04 


9. 


3.709 


4.622 


1.794E+47 


4.303E+46 


1.301E+41 


3.785E+47 


0.240E+00 


0.725E-06 


5. 


2.870 


4.440 


1.097E+45 


8.629E+41 


7.605E+36 


3.760E+46 


0.787E-03 


0.693E-08 



The ionising properties of ZAMS models for masses 
between 5 and 1000 Mq are given in Table ||^. 

The following differences are obtained between vari- 
ous atmosphere models for T^s between 50 and 100 kK. 
The plane parallel TLUSTY models show some variation 
of g(He+) at 70 kK with gravity, with low logg leading to 
higher q(He+). As expected, for very hot stars the CMF- 
GEN photosphere-wind models (filled triangles in Fig. ||) 
show little deviation from the plane parallel models, since 
the He ii continuum is formed in progressively deeper lay- 
ers of the atmosphere (e.g. Husfeld et al. 1984, Clegg & 
Middlemass 1987). The two models at Toff ^ 70 kK (cor- 
responding to a ZAMS mass of 25 M0 ) show an increase 
of g(He+) of up to - 1 dex with respect to the TLUSTY 
model of same logg, due to a depopulation of the He 11 
groundstate induced by the stellar wind (Gabler et al. 

® The luminosities and T^ff are taken from the set of tracks 
without mass loss, complemented by the 1000 Mq model of 
Klapp (1983). As expected the differences between the vari- 
ous ZAMS models is small. The approximate stellar models of 
Bromm et al. (2001) and TSOO yield somewhat higher Teff by 
up to ~ 6 kK and 30 kK respectively. 



1989). Towards lower Toff differences between static and 
spherically expanding atmospheres becomes progressively 
more important for (/(He^). However, in the integrated 
populations considered later, the emission of He 11 ionis- 
ing photons will be dominated by the hottest objects. We 
are therefore likely left with the situations where either 
possible wind effects have no significant impact on the 
total q{lle^)/q(H) hardness, or q{lle'^)/q(H) is too small 
to have observable consequences. In conclusion, it appears 
safe to rely only on non-LTE plane parallel atmospheres to 
properly describe the ionising fluxes of metal-free stellar 
populations. 

Can other effects, not included in the present model at- 
mosphere grids, affect the ionising spectra ? The effect of 
non-coherent electron scattering could be important for 
Toff ^ 40 kK (Rybicki & Hummer 1994). Test calcula- 
tions with the CMFGEN code of Hilher & Miller (1998) 
for some of the above CMFGEN models indicate negligi- 
ble changes on the ionising photon fluxes. Compton scat- 
tering, instead of the commonly implemented Thomson 
scattering, should not affect the ionising fluxes of normal 
Pop HI stars, as differences appear only above Toff 150 
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Table 4. Time averaged properties calculated from Pop III tracks with no mass loss 

Mini lifetime Q(H) Q(He°) Q(He+) ^(Ha) Q(He°)/Q(H) Q(He+)/Q(H) 



1000. not available 



500 


00 


1 


899E+06 


6 


802E+50 


3 


858E+50 


5 


793E+49 


7 


811E+50 


0.567E+00 


0.852E-01 


400 


00 


1 


974E+06 


5 


247E+50 


3 


260E+50 


5 


567E+49 


5 


865E+50 


0.621E+00 


0.106E+00 


300 


00 


2 


047E+06 


3 


754E+50 


2 


372E+50 


4 


190E+49 


4 


182E+50 


0.632E+00 


0.112E+00 


200 


00 


2 


204E+06 


2 


624E+50 


1 


628E+50 


1 


487E+49 


2 


918E+50 


0.621E+00 


0.567E-01 


120 


00 


2 


521E+06 


1 


391E+50 


7 


772E+49 


5 


009E+48 


1 


608E+50 


0.559E+00 


0.360E-01 


80 


00 


3 


012E+06 


7 


730E+49 


4 


317E+49 


1 


741E+48 


8 


889E+49 


0.558E+00 


0.225E-01 


60 


00 


3 


464E+06 


4 


795E+49 


2 


617E+49 


5 


136E+47 


5 


570E+49 


0.546E+00 


0.107E-01 


40 


00 


3 


864E+06 


2 


469E+49 


1 


316E+49 


8 


798E+46 


2 


903E+49 


0.533E+00 


0.356E-02 


25 


00 


6 


459E+06 


7 


583E+48 


3 


779E+48 


3 


643E+44 


9 


387E+48 


0.498E+00 


0.480E-04 


15 


00 


1 


040E+07 


1 


861E+48 


8 


289E+47 


1 


527E+43 


2 


526E+48 


0.445E+00 


0.820E-05 


9 


00 


2 


022E+07 


2 


807E+47 


7 


662E+46 


3 


550E+41 


5 


576E+47 


0.273E+00 


0.126E-05 


5 


00 


6 


190E+07 


1 


848E+45 


1 


461E+42 


1 


270E+37 


6 


281E+46 


0.791E-03 


0.687E-08 



Table 5. Time averaged properties calculated from Pop III tracks with strong mass loss 



Mini 


lifetime 


Q(H) 


(3(He°) 


Q(He+) 


0(H2) 


Q(He°)/Q(H) 


g(He+)/0(H) 


1000. 


2.430E+06 


1.863E+51 


1.342E+51 


3.896E+50 


2.013E+51 


0.721E+00 


0.209E+00 


500. 


2.450E+06 


7.719E+50 


5.431E+50 


1.433E+50 


8.345E+50 


0.704E+00 


0.186E+00 


300. 


2.152E+06 


4.299E+50 


3.002E+50 


7.679E+49 


4.766E+50 


0.698E+00 


0.179E+00 


220. 


2.624E+06 


2.835E+50 


1.961E+50 


4.755E+49 


3.138E+50 


0.692E+00 


0.168E+00 


200. 


2.628E+06 


2.745E+50 


1.788E+50 


2.766E+49 


3.028E+50 


0.651E+00 


O.lOlE+00 


150. 


2.947E+06 


1.747E+50 


1.156E+50 


2.066E+49 


1.917E+50 


0.662E+00 


0.118E+00 


100. 


3.392E+06 


9.398E+49 


6.118E+49 


9.434E+48 


1.036E+50 


0.651E+00 


O.lOOE+00 


80. 


3.722E+06 


6.673E+49 


4.155E+49 


4.095E+48 


7.466E+49 


0.623E+00 


0.614E-01 



kK (Hubeny et al. 2001). X-rays, originating from stellar 
wind instabilities or interactions with magnetic fields, can 
potentially increase the high energy (He"*") ionising flux in 
stars with weak winds (MacFarlane et al. 1994). For Pop 
III stars such a hypothesis remains highly speculative. The 
adopted model atmospheres should thus well describe the 
spectra of the objects considered here. 



3.1.1. Comparisons with other Pop III calculations 

As expected, our calculations for individual stars of a given 
Toff agree quite well with those of TSOO based also on 
the TLUSTY model atmospheres. As briefly mentioned 
above, small differences are found with their ZAMS po- 
sitions calculated with simplified stellar models. For in- 
tegrated populations the use of a proper ZAMS leads to 
somewhat softer spectra (cf. Sect. 4.1). 

As pointed out by Bromm et al. (2001) the spectral 
properties of stars with M ^ 300 Mq are essentially in- 
dependent of stellar mass, when normalised to unit stellar 
mass. The degree to which this holds can be verified from 
Table ^. If correct, it would imply that the total spectrum 
of such a population would be independent of the IMF and 
only depend on its total mass. Our results for the 1000 Mq 



ZAMS star are in rather good agreement with the ionising 
fluxes of Bromm et al. (2001, their Table 1) ^. However, 
as already apparent from their Fig. 3 the deviations from 
this approximate behaviour are not negligible for the He^ 
ionising flux. In consequence we find e.g. Q(II) reduced by 
12 % and He ii recombination line strengths reduced by 
a factor 2 for a population with a Salpeter IMF from 300 
to 1000 Mq! Larger difference are obviously obtained for 
IMFs with Miow < 300 Mq. Further differences (due to the 
inclusion of nebular continuous emission) with Bromm et 
al. (2001) are discussed in Sect. |[ 



3.2. Time averaged ionising properties 

As clear from Fig. |l|, massive stars evolve rapidly toward 
cooler temperatures during their (short) lifetime in the 
absence of strong mass loss. It is thus evident that espe- 
cially the He ii ionising flux will strongly decrease with 
age. To quantify this effect we have calculated the lifetime 



^ Our values for Q(H), Q(He°), and Q(He+) differ by 0, +3, 
and - 25 % respectively from the values of Bromm et al.. 
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Fig. 5. Spectral energy distribution (SED) including H and He recombination lines for model B (solid line, see Table 
|l| for line identifications). Left panel: ZAMS population. The pure stellar continuum (neglecting nebular emission) 
is shown by the dashed line. For comparison the SED of the Z = 1/50^0 population (model ZL: Salpcter IMF from 
1 - 150 Mq) is shown by the dotted line. The vertical dashed lines indicate the ionisation potentials of H, IIe°, and 
IIe+. Note the presence of the unique He ii features (shown as thick dashed lines) and the importance of nebular 
continuous emission. Right panel: temporal evolution of the spectrum for 0, 2 and 4 Myr showing the rapid change 
of the emission line spectrum, characterised by the disappearance of the He ii lines 



averaged ionising flux Qi along the evolutionary tracks, 
defined by 

r«*(M) 



.(M) = 



Jo 



^it,M)dt 



(5) 



U{M) 

where is the stellar lifetime. The results calculated for 
both sets of tracks (no mass loss, strong mass loss) are 
given in Tables || and || respectively. For convenience we 
also provide simple fits (linear or quadratic) of these values 
as a function of the initial mass M in Table |^. Together 
with ti, the average ionising flux Qi defined in this manner 
in particular allows one to calculate the correct values of 
the ionising flux, recombination line luminosities etc. for 
integrated populations at equilibrium, i.e. for the asymp- 
totic value obtained for a constant star formation rate. 

The time averaged hardness as a function of mass is 
compared to the ZAMS value in Fig. |3[ For tracks without 
mass loss, Q(He+)/(3(H) is found to be a factor of two 
or more below the ZAMS value. Smaller differences are 
obtained for obvious reasons for Q(Re^)/Q(R). Only in 
the case of very strong mass loss, the equilibrium value 
of the hardness of the ionising spectrum is found to be 
essentially identical as the ZAMS value. 

For Pop I stars the effect of time dependent and time 
averaged ionising fluxes and the spectral hardness based 
on recent atmosphere models have been studied earlier in 
various contexts (analysis of O star populations from in- 



tegrated spectra, studies of the diffuse ionised gas; see e.g. 
Schaerer 1996, 1998). For comparison sake the resulting 
fits for solar metallicity and 1/50 Zq are also given in 
Table |. 

4. Population III "galaxies" 

4.I. Time evolution of integrated spectra 

The temporal evolution of the integrated ionising photon 
flux of an instantaneous burst and the case of constant 
star formation is plotted in Fig. ^ for different IMFs. Most 
notable is the rapid decrease of the He ii ionising flux in 
a burst, due to the evolution of the most massive stars off 
the hot ZAMS and due to their short lifetime. 

Our values compare as follows with other calculations 
of Q(H) available in the literature: when scaled to the 
same lower mass cut-off of the IMF (Miow = 0.024 Mq), 
our ZAMS model for IMF A (Table |) agrees with Q(H) 
(age ~ 0) of Cojazzi et al. (2000). The same holds for (5(H) 
of Ciardi et al. (2000), when adopting their IMF (Salpeter 
from 1 to 40 Mq). A good agreement is found with (5(H) 
of TSOO (ZAMS population with a Salpeter IMF from 0.1 
to 10 Mq), whereas the slightly cooler ZAMS (cf. Sect. 
p.l.lD implies a somewhat softer spectrum ((5(He"'")/(5(H) 
~ 0.03 instead of 0.05). Differences with Bro mm et al. 
(2001) have already been discussed above (Sect. 3.1.1 ). 
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0.6 



0.4 



0.3 



q(He)/q(H) 



ZAMS 




Fig. 3. Hardness of the ionising spectrum expressed by 
g(He°)/(7(H) and g(He+)/q(H) as a function of the initial 
stellar mass. Dashed lines show the hardness of all ZAMS 
models. The thick lines (dashed: no mass loss, dotted: high 
mass loss tracks) show the average spectral hardness when 
stellar evolution effects are taken into account. Note in 
particular the importance of stellar evolution leading to a 
decrease of the hardness of the He ii ionising spectrum by 
approximately a factor of 2 

Spectral energy distributions of integrated zero metal- 
licity stellar populations are shown in Fig. ||for the case of 
a Salpeter IMF from 1 to 500 Mq and instantaneous bursts 
of ages (ZAMS), 2, and 4 Myr. Overplotted on the con- 
tinuum (including stellar + nebular emission: solid lines) 
are the strongest emission lines for illustration purpose. 
In the left panel we show for comparison the spectrum of 
a burst at low metallicity (1/50 Zq, Salpeter IMF from 
1 - 150 M0; dotted line). The striking differences, most 
importantly in the ionising flux above the He 11 edge (> 54 
eV), have already been discussed by TSOO. The compari- 
son of the total spectrum (solid line) with the pure stel- 
lar emission (dashed) illustrates the importance of nebu- 
lar continuous emission neglected in earlier studies (TSOO, 
Bromm et al. 2001), which dominates the ZAMS spectrum 
at A > 1400 A. The nebular contribution, whose emission 
is proportional to Q(H) (Eq. ^, depends rather strongly 
on the age, IMF, and star formation history. For the pa- 
rameter space explored here (cf. Table we find that 
nebular continuous emission is not negligible for bursts 
with ages ^ 2 Myr and for constant star formation mod- 
els. 

The spectra in Fig. ^ show in addition to the H and 
He I recombination lines the presence of the strong He 11 



Table 6. Least square fits to the time averaged properties 
for Pop III stars, metal-poor objects {Z = 1/50 Zq) and 
solar metallicity stars for comparison. The time averaged 
quantities \ogy{M) — uq + aix + a2X^ + a^x^ (col. 1) are 
fitted as a function of the initial mass {x = log{M / Mq)) as 
2nd or 3rd order polynomials over the mass interval given 
in col. 2. For masses lower than the range considered the 
corresponding photon flux Q is negligible. Q values are in 
units of photon s^^, stellar lifetimes t^, in yr. 



Quantity y 


mass range 


ao 


ai 


02 


as 


Z — tracks with no mass loss: 








Q(H) 


9 - 500 Mq 


43.61 


4.90 


-0.83 




Q(H) 


5 - 9 M© 


39.29 


8.55 






Q(He°) 


9 - 500 Mo 


42.51 


5.69 


-1.01 




Q(He«) 


5 - 9 M0 


29.24 


18.49 






Q(He+) 


5 - 500 Mq 


26.71 


18.14 


-3.58 




Q(H2) 


5 - 500 Mo 


44.03 


4.59 


-0.77 




i* 


5 - 500 Mo 


9.785 


-3.759 


1.413 


-0.186 


Z — tracks with strong mass 


loss: 








Q(H) 


80 - 1000 Mq 


46.21 


2.29 


-0.20 




Q(He«) 


80 - 1000 Mq 


45.71 


2.51 


-0.24 




Q(He+) 


80 - 1000 Mq 


41.73 


4.86 


-0.64 




Q(H2) 


80 - 1000 Mq 


46.25 


2.30 


-0.21 






80 - 1000 Mq 


8.795 


-1.797 


0.332 




Solar metallicity tracks: 










Q(H) 


7 - 120 Mo 


27.89 


27.75 


-11.87 


1.73 


Q(He«) 


7 - 120 Mo 


1.31 


64.60 


-28.85 


4.38 


Q(H2) 


15 - 120 Mq 


41.90 


7.10 


-1.57 




U 


7 - 120 Mo 


9.986 


-3.497 


0.894 




Z = 1/50 Zq 


tracks: 










Q(H) 


7 - 150 Mo 


27.80 


30.68 


-14.80 


2.50 


g(He°) 


20 - 150 Mq 


16.05 


48.87 


-24.70 


4.29 


g(He+) 


20 - 150 Mq 


34.65 


8.99 


-1.40 




Q(H2) 


12 - 150 M© 


43.06 


5.67 


-1.08 




U 


7 - 150 Mo 


9.59 


-2.79 


0.63 





AA 1640, 3203, and 4686 recombination lines, which — 
due to the exceptional hardness of the ionising spectrum 
— represent a unique feature of Pop III starbursts com- 
pared to metal enriched populations (cf. TGSOl, Oh ct al. 
2001a, Bromm et al. 2001). Another effect highhghted by 
this Figure is the rapid temporal evolution of the recom- 
bination line spectrum. Indeed, already <; 3 Myr after the 
burst, the high excitation lines are absent, for the reasons 
discussed before. In the case of constant star formation, 
the emission line spectrum at equilibrium is similar to a 
burst of ~ 1.3 Myr (cf. Fig. |). 

Figure 1^ shows the temporal evolution of the line inten- 
sity of He II A1640 with respect to H/3 for different IMFs, 
instantaneous bursts (upper panel) or constant star for- 
mation (lower panel), and the two sets of stellar tracks 
(see model designations in Table ||) . As the lifetime of the 
stars affected by the IMF variations considered here is ~ 
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Fig. 4. Ionising photon fluxes (9(H), Q(He^) and Q(He+) 
as a function of time for instantaneous burst models (up- 
per panel) or constant star formation (lower panel). For 
each group of 3 lines (using the same symbols) Q(H), 
Q(IIe'') and (9(IIe+) are given from top to bottom. Dif- 
ferent IMFs are illustrated by the dotted (model ID A), 
solid (model B) and long-dashed lines (model C) respec- 
tively (cf. Table H). The units are photon s~^ M©"^ for 
the burst and photon s~^ normalised per unit SFR (in 
Mq yr~^) for constant SF. Note in particular the rapid 
decrease of the He ii ionising flux ~ 3-4 Myr after the 
burst 



instantaneous burst 




2 4 6 8 10 

age [Myr] 



Fig. 6. Temporal evolution of the relative He ii A1640/H/3 
line intensity for various models. Upper panel: instanta- 
nous burst models with different IMFs and for no mass loss 
tracks (dotted, dashed, sohd) and high mass loss (dash- 
dotted, long dash ~ dotted; see symbols in Table ^). Lower 
panel: same as above for constant star formation models. 
The shaded ranges on the right indicate He ii A1640/H/? 
intensities of quasars and maximum values expected for 
metal-poor H ii regions. For the IMFs considered here 
/(1640) //(H/3) can reach the high values typical of quasars 
for very young bursts or in the (extreme) case of high mass 
loss. 



2-4 Myr, no changes are seen at ages 4 Myr. The 
following points are worthwhile noticing: 

1 ) aheady for a "normal" IMF with M^p = 100 M© the 
maximum He ii A 1640 intensity is exceptionally large 
for pure stellar photoionisation compared to metal- 
poor stellar populations (cf. TSOO); 

2) extending the IMF up to 500 - 1000 Mg leads to 
/(1640)//(H/3) up to - 1.5 - 2; 

3) in the unlikely case of high mass loss the He ii A1640 
intensity is maintained for up to 4 Myr; 

4) for models without mass loss, the equilibrium val- 
ues of /(1640)//(H/3) are of the order of ~ 0.1 - 
0.6 for the IMF range considered here. For compari- 
son, /(1640)//(H^) - 1.0 - 3.2 for a power law spec- 
trum Ln oc and typical values a ~ 1 - 1.8 for 
quasars (Elvis et al. 1994, Zheng et al. 1997), and the 
largest intensities in (rare) metal-poor H ii regions are 
/(1640)//(H/3) <, 0.15, as expected from their mea- 
sured He II A4686 intensity. 



The equivalent widths for the strongest H and He ii 
lines and their dependence on the IMF and stellar tracks 
are shown in Figs. ^ and |^. Maximum equivalent widths of 
~ 700 and 3200 A are predicted for the optical lines of H/3 
and Ha (not shown) respectively. W(Lyman-a) reaches 
values up to ~ 2100 A at the ZAMS, but is potentially 
affected or destroyed by radiation transfer or other effects 
(cf. Tenorio-Tagle et al. 1999). The maximum equivalent 
width of the He ii AA 1640 and 4686 recombination lines on 
the ZAMS is predicted to be ~ 100 and 150 A respectively. 

What variations are typically obtained for different 
IMF slopes ? Varying a between 1. and 3. one finds the 
following changes for bursts at age Myr with respect 
to the Salpeter slope: for models A and B /(1640)//(H/3) 
(and the other He ii intensities) varies between ^ - 60 % 
and + 50 %, W^(Ly— a) changes by not more than ~ 20 
%, and W^(He ii A1640) varies between ~ - 50 % and + 60 
%. For obvious reasons, populations most biased towards 
massive stars (i.e. large values of Miow) are the least sensi- 
tive to the exact IMF slope (cf. Bromm et al. 2001). E.g., 
the above quantities vary by 20 % for model C. 
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2500 




200 



Fig. 7. Temporal evolution of Lyman-a (top) and H/3 
(bottom) for instantaneous burst models. Note that very 
large maximum equivalent widths are obtained. As the 
IMF and track variations considered here concern only the 
most massive stars (M 50 - 100 M0) the predictions 
are only affected at young ages (^ 4 Myr) 



Much larger Lyman-a and He II A1640 equivalent 
widths (^ 3100 and 1100 A respectively) have been pre- 
dicted for a ZAMS population of exclusively massive stars 
by Bromm et al. (2001). This overestimate is due to two 
effects acting in the same direction: their simplified use of 
the spectrum of a 1000 Mq star with the hardest spec- 
trum as representative for the entire population^, and 
most importantly the neglect of nebular continuous emis- 
sion, which — for the case of a Salpeter IMF from 300 
to 1000 Mq — contributes ~ 75 % of the total light at 
A 1640 and thus strongly reduces the He 11 A1640 equiva- 
lent width (cf. above). Accounting for these effects yields 
W(Lyman-a) - 2000 and W(1640) - 120 A for their IMF. 

Our calculation of the nebular continuous emission as- 
sumes density bounded objects, i.e. negligible escape of 
Lyman continuum photons out of the "Pop HI galaxies" 
(/esc ^ 1). If this is not the case, one likely has a situation 
where nevertheless all He 11 ionising photons are absorbed 
— i.e. L(He 11 A1640) remains identical — while some H 
ionising photons escape thereby reducing the continuum 
emission (proportional to (3(H) (Eq. Although the es- 
cape fraction of ionising photons in distant galaxies or the 
first building blocks thereof remains badly known, there 
are theoretical and observational indications for rather low 




For a Salpeter IMF from 300 - 1000 Mq this leads to an 
overestimate of L(1640) by a factor of 2 (cf. Sect. H). 



Fig. 8. Same as | for the He 11 A1640 (top) and He 11 A4686 
(bottom) recombination lines. Note that large maximum 
equivalent widths. 



escape fractions, of the order of ^ 10 % (see e.g. Leitherer 
et al. 1995, Dove et al. 2000, Steidel et al. 2001, Hui et 
al.2001). For such moderate an escape fraction, nebular 
continuous emission remains thus dominant when strong 
emission lines are expected, leading to He 11 equivalent 
widths not much larger than predicted here. Estimates 
for small Pop HI halos suggest, however, large escape frac- 
tions (Oh et al. 2001b). In short, the caveat regarding the 
uncertainty on /esc, which affects both emission line lu- 
minosities and nebular continuum emission, must be kept 
in mind. A more detailed treatment including geometrical 
effects, radiation transfer etc. (cf. Ciardi et al. 2001) are 
beyond the scope of our estimates. 

Regarding the equivalent widths, we conclude that 
in any case the maximum values of W^(Lyman-a) and 
PF(He II A1640) of metal-free populations are ~ 3-4 and ~ 
80 times larger than the values predicted for very metal- 
poor stellar populations {Z = 1/50 Zq) with a "normal" 
IMF including stars up ~ 100 Mq. When measured in 
the rest-frame, these will further be amplified by a factor 
{1 + z), where z is the redshift of the source. 

4-2. Star formation indicators 

In the case of constant star formation, the ionising prop- 
erties tend rapidly (typically over ~ 5 - 10 Myr depending 
on the IMF) to their equilibrium value, which is then pro- 
portional to the star formation rate (SFR). In particular. 
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recombination line luminosities L/ are then simply given 

by 

i--«(i-/o.)<3.(j|^)=/,(-i^). 

The predicted values /; (in erg s~^) for He il A1640 and 
Ha for the IMFs considered here are given in Table |^ for 
the case of a negligible escape fraction of ionising photons 
(/esc = 0). For other lines, or to compute the Qi not listed 
explicitly, the corresponding value can easily be derived 
using Eq. ^ and the emission coefficients in Table |^. Also 
listed are the He i ionising production and the photon flux 
in the Lyman Werner band of H2 for a star formation rate 
of 1 M0 yr^^. For comparison with SFR indicators used 
for "normal" galaxies we also indicate the mass conver- 
sion factor Cm expressing the relative masses between our 
model IMF and a "standard" Salpeter IMF with Miow = 
0.1 and Mup = 100 M© used frequently throughout the 



literature, i.e. cm = J^""^ M^{M)dm/ J^^^^^ M^{M)dm. 

Compared to the "standard" Ha SFR indicator for so- 
lar metallicity objects (cf. Kcnnicutt 1998, Schaerer 1999) 
L(Ha) is ~ 2.4 times larger for a Pop HI object with the 
same IMF, given the increased Q(H) production (cf. TSOO 
and above). Obviously, for IMFs more weighted towards 
high mass stars, the SFR derived from Ha (or other H 
recombination lines) is even more reduced compared to 
using the standard Pop I SFR indicator. If detected. He 11 
A1640 can also be used as a star formation indicator as al- 
ready pointed out by TGSOl. Their corresponding value is 
/i640 = 8.4 1 0"*° X /cvoi, with /cvoi ~ 0.4 (2.) accounting in 
an approximate way for stellar evolution effects without 
(with) mass loss. Including the full sets of evolutionary 
tracks we find /i64o = 1-91 10"*", for the 1-100 Mq IMF, 
i.e. less He 11 emission than TGSOl. Varying the IMF slope 
a between 1. and 3. leads to changes of /hq (/i64o) by ~ ± 
0.5 (0.9) dex for models A and B, and to minor changes 
(^ 0.1 dex) for model C. With the fits given in Table || 
the line luminosities can easily be computed for arbitrary 
IMFs. 

In passing we note that the He"*" ionising flux per unit 
SFR used by Oh et al. (2001a) to estimate the expected 
number of sources with detectable He 11 recombination 
lines has been overestimated: for a constant SFR and a 
Salpeter IMF up to 100 M© the hardness (5(He+)/Q(H) 
~ 0.005 is a factor 10 lower than their Q value adoptecj^. 
However, the Lyman continuum production is increased 
over the value they adopted, leading thus to a net reduc- 
tion of their flux Jf^^^ by a factor 1.9 (4.8) for a Salpeter 
IMF with Miow = 1 (0.1) and Mup=100 M© (cf. Tables 
1^ and 0) . Further issues regarding the photometric prop- 
erties and the detectability of Pop HI sources will be dis- 
cussed in Pello & Schaerer (2001). 



J^.S. Photon and heavy element production 

Various other properties related to the photon and heavy 
element production are given in Table | for the case of 
constant star formation normalised to SFR = 1 M© yr~^. 
Columns 5 and 6 give (5(He°), the He" ionising flux, and 
the photon flux in the Lyman- Werner bands Q(LW) re- 
spectively, both in units of (photon s~^)/(Mq yr^^). The 
average energy Eq (in eV) of the Lyman continuum pho- 
tons is in col. 7. 

To estimate the metal production of Pop HI objects we 
have used the SN yields of Woosley & Weaver (1995) for 
progenitor stars between 8 and 40 M© and the results of 
Hegcr ct al. (2000) and Heger & Woosley (2001) for very 
massive objects of pair creation SN originating from stars 
in the mass range ~ 130 to 260 M© (cf. also Ober et al. 
1993). As the pair creation SN models of Heger and col- 
laborators are calculated starting with pure helium cores, 
we have adopted the relation Ma ~ 0.5 Mini from Ober 
et al. (1983) to calculate the initial mass Mini. We neglect 
the contribution from longer lived intermediate mass stars 
(1 ^ M /Mq ^ 8). Non-rotating stars more massive than 
~ 260 M© are expected to collapse directly to a black hole 
producing no metals (cf. Rackavy et al. 1967, Bond et al. 
1984). 

The SN rate SNR including "normal" SN and pair cre- 
ation SN is given in col. 8 of Table |^ (in units of SN per 
solar mass of stars formed). The total mass of metals Mej 
produced including all elements heavier than He is given 
in col. 9. Columns 10, 11, and 12 indicate the contribution 
of C, O and Si. The ejected masses are given in units solar 
mass per unit SFR. Note that the models C and E with 
a lower mass cut-off Miow = 50 M© of the IMF represent 
cases with solely pair creation SN (no SN type II) , models 
A (plus the cases ZS and ZL with non-zero metallicity) 
include only SNII, and models B and D include both SN 
types. Finally the dimensionless "conversion efficiency" of 
ionising photons to rest mass (cf. Madau & ShuU 1996) 

„ _ /o°° C ^-'^'"^^ _ Q(H) * Eq 



M,^c? 



M.ic-' 



® Oh et al. (2001a) adopt the ZAMS population of TSOO 
as representative of a stellar population at equilibrium with a 
given SFR. 



with the spectral luminosity L^, is given in col. 13. It is 
useful to remind that, since both the metal and photon 
production is dominated by massive stars, 77 is indepen- 
dent of the IMF at masses M ^ 5-8 M©, and depends 
little on the exact shape of the IMF. 

The SNR per unit mass depends little on the upper 
mass cut-off, as most of the mass resides in low mass stars. 
For cases C and E with Miow = 50 M© the SNR is reduced, 
since pair creation SN from stars of a limited mass range 
(~ 130 - 260 M©) are the only explosive events. 

As evident from the yields of Woosley & Weaver (1995) 
and shown in Table |7[ Population HI stars convert - for 
"standard" IMFs - less of their initial mass into metals 
than stars of non-zero metallicity^. However, for IMFs 

® Mej and rj given in Table ^ for solar metallicity (model ZS) 
have been calculated with the code of Cervino et al. (2000) 
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Table 7. Line luminosities, supernova rates, and photon 
normalised to SFR = 1 Mq yr~^ 

ID CM /Ha /l640 Q(He°) Q(LW) 

[erg s-i] [(ph s-1)/(Mq yr"!)] 



and metal production for constant star formation models 

Eq SNR Mej Moj(C) Moj(O) Moj(Si) »7 
[eV] [M0-1] [Mq/(Mq yr-i)] 



A 2.55 
B 2.30 
C 14.5 



7.88e+41 
1.03e+42 
3.32e+42 



1.91e+40 
9.98e+40 
8.38e+41 



3.15e+53 
4.34e+53 
1.54e+54 



3.17e+53 
3.29e+53 
4.52e+53 



26.70 
27.86 
29.69 



0.019 
0.016 
0.0016 



0.0076 

0.030 

0.14 



0.0013 
0.0023 
0.0065 



0.0030 

0.014 

0.067 



0.0005 0.065 
0.0045 0.022 
0.025 0.016 



D 2.30 1.14e+42 3.12e+41 5.26e+53 3.36e+53 30.55 0.016 
E 12.35 4.10e+42 2.33e+42 2.21e+54 3.65e+53 35.15 0.0013 



Solar metallicity tracks: 

ZS 2.55 3.16e+41 5.12e+39 5.56e+52 3.60e+53 20.84 0.019 0.038* 
Z = 1/50 Zq tracks: 

ZL 2.47 6.92e+41 1.70e+39 1.556+53 5.31e+53 21.95 0.019 0.023 



0.0025 0.014 



0.0036* 
0.0014 0.014 



Including stellar wind mass loss, SN Ibc, and SNII 



favouring more strongly massive stars, allowing in partic- 
ular for pair creation SN ejecting up to half of their initial 
mass, the metal production per unit stellar mass can be 
similar to or larger than for solar metallicity. 

The photon conversion efficiency of rest mass to ion- 
ising continuum, 77, of metal-poor and metal-free popula- 
tions is increased by factors of ^ 3 - 18 compared to solar 
metallicity. In great part this is due to increased ionising 
photon production. The underlying changes of all quanti- 
ties on which rj depends are listed in Table 0. The increase 
of 77 from Zq to 1/50 Z©]^ indicates already a substan- 
tially larger photon to heavy element production which 
must have occurred in the early Universe. 

In passing we note that our value of 77 '--^ 0.36% (or 0.5 
% neglecting stellar mass loss and SNIbc) is larger than 
the value of 0.2 % given by Madau & ShuU (1996), which, 
given the quoted assumptions appears to be erroneous. 
This can be easily seen by noting the good agreement 
between our photon production calculation and other cal- 
culations in the literature (e.g. the Ha SFR calibration 
of Kennicutt 1998, Glazebrook et al. 1999). In addition 
one can easily estimate the metal production from noting 
that Moj(M) = /Mini, with typically / 0.1 (Woosley & 
Weaver 1995) and integrating over the IMF, which con- 
firms our calculations. 

Finally the predicted ejecta of carbon, oxygen, and sil- 
icon show some interesting variations with the IMF. From 
the "standard" IMF (A) to an IMF producing exclusively 
very massive stars (model C), the oxygen/carbon ratio 
(0/C) ~ 0.8 (0/0)0 increases by a factor ~ 4, while in- 
taking into account stellar wind mass loss, SN Ibc, and SNII. 
Neglecting mass loss and assuming only type II SN we obtain 
Mej = 0.027 Mq and 77 = 0.005, differing by ~ 40 % from the 
tabulated value. 

Assuming Mup = 100 Mq as for models A and ZS reduces 
Q(H) by 16 %. 



creases from Si/0 ^ 1.8 (Si/O)© by a factor of ^ 10! It 
is interesting to compare these values with current mea- 
surements of abundance ratios in the Lyman-a forest, 
which indicate an overabundance of Si/0 2 - 3 (Si/O)© 
(Songaila & Oowie, 1996, Oiroux & ShuU 1997), but possi- 
bly up to ~ 10 times solar, depending on the adopted ion- 
ising UV background radiation field (Savaglio et al. 1997, 
Giroux & ShuU 1997). If the ejecta from Pop III objects are 
responsible for the bulk of metals in the Lyman-a forest, it 
would appear that IMFs strongly favouring massive stars 
seem to be excluded by studies finding a modest Si/O over- 
abundance. The use of more detailed chemical evolution 
models including also intermediate mass stars (cf. Abia 
et al. 2001) and additional studies on abundances in the 
Lya forest should hopefully provide stronger constraints 
on early nucleosynthesis and the IMF of Pop III stars. 

5. Summary and discussion 

The main aim of present study was to construct realistic 
models for massive Population III stars and stellar popu- 
lations to study their spectral properties, including their 
dependence on age, star formation history, and most im- 
portantly on the IMF, which remains very uncertain for 
such objects. 

We have calculated extensive sets of non-LTE model 
atmospheres appropriate for metal-free stars with the 
TL US TY code (Hubeny & Lanz 1995) for plane paral- 
lel atmospheres. The comparison with non-LTE models 
including stellar winds, constructed with the comoving 
frame code CMFGEN of HiUier & Miller (1998), shows 
that even in the presence of some putative (weak) mass 
loss, the ionising spectra of Pop III stars differ negligi- 
bly from those of plane parallel models, in stark contrast 
with Pop I stars (Gabler et al. 1989, Schaerer & de Koter 
1997). As aheady discussed by Tumhnson & ShuU (2000), 
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the main salient property of Pop III stars is their increased 
ionising flux, especially in the He"*" continuum (> 54 eV). 

The model atmospheres have been introduced together 
with recent metal free stellar evolution tracks from the 
Geneva and Padova groups (Feijoo 1999, Desjacques 2000, 
Marigo et al. 2000) and older Pop III tracks assuming 
strong mass loss (Klapp 1983, El Eid et al. 1983) in the 
evolutionary synthesis code of Schaerer & Vacca (1998) to 
study the temporal evolution of individual Pop III stars 
and stellar populations. 

The main results obtained for individual Pop III stars 
are the following (Sect. ||): 

— Due to their redward evolution off the zero age main se- 
quence (ZAMS) the spectral hardness measured by the 
He"*"/!! ionising flux is decreased by a factor ^ 2 when 
averaged over their lifetime. If such stars would suffer 
strong mass loss, their spectral appearance could, how- 
ever, remain similar to that of their ZAMS position. 

— The ionising photon fluxes of H, He°, He"*", and the flux 
in the Lyman- Werner band dissociating H2 have been 
tabulated and detailed fit formulae are given for stars 
from 5 to 1000 M0. These can in particular be used to 
calculate with arbitrary IMFs the properties of ZAMS 
populations or cases of constant SFR. For comparison 
the same quantities are also given for solar metallicity 
and 1/50 Zq. 

The integrated spectral properties including line emis- 
sion from H, He i. He 11 recombination lines (treated e.g. 
by TSOO, Bromm et al. 2001) and nebular continuous emis- 
sion — neglected in all earlier studies — have been cal- 
culated for instantaneous bursts and constant formation 
rates. Various assumptions on the mass cut-offs of the IMF 
including a "standard" IMF from 1 to 100 Mq or IMFs 
forming exclusively very massive stars have been consid- 
ered (see Table ^). 

The main results regarding integrated stellar popula- 
tions are as follows (Sect. ^): 

— For young bursts and the case of a constant SFR, neb- 
ular continuous emission dominates the spectrum or 
is at least important to properly predict it redward 
of Lyman-a if the escape fraction /esc of ionising pho- 
tons out of the observed region is small or negligible. In 
consequence predicted emission line equivalent widths 
are considerably smaller than found in earlier stud- 
ies, whereas the detection of the continuum is eased. 
Maximum equivalent widths of W(Lya) ~ 1700 and 
W(He II A1640) ~ 120 A are predicted for a Salpeter 
IMF with Mup - 500 - 1000 Mq. Nebular line and con- 
tinuous emission strongly affect the broad band photo- 
metric properties of Pop III objects and (see Schaerer 
& Pello 2001, Pello & Schaerer 2001). Detailed treat- 
ment models including geometrical effects, radiation 
transfer etc. (cf. e.g. Ciardi et al. 2001) will be useful 
to constrain /osc more precisely. 



— Due to the redward stellar evolution and short lifetimes 
of the most massive stars, the hardness of the ionising 
spectrum decreased rapidly, leading to the disappear- 
ance of the characteristic He 11 recombination lines af- 
ter ~ 3 Myr in instantaneous bursts. 

— He II A1640, Ha (and other) line luminosities usable as 
indicators of the star formation rate are given for the 
case of a constant SFR (Table ^). For obvious reasons 
such indicators depend strongly on the IMF. 

— Accounting for SNII and pair creation supernovae we 
have also calculated supernova rates, the metal pro- 
duction, and the conversion efficiency of ionising pho- 
tons to rest mass, 77 (Table 0). For a "standard" IMF 
with Mup = 100 M0 the increased photon production 
and the reduced metal yields of Pop III stars lead to 
an increase of 77 by a factor of 5 to 18 when compared 
to low (1/50 Zq) or solar metallicity. In the presence 
of very massive stars (M ~ 130 - 260 Mq) leading to 
pair creation SN, somewhat lower values of ry ~ 1.6 to 
2.2 % are obtained. In any case metal-poor and metal- 
free populations are much more efficient producers of 
ionising photons "per baryon" (per unit heavy element 
rest mass energy) than at solar metallicity, where rj ~ 
0.36 %. 

— Finally the ejected masses of C, O, and Si of Pop III 
objects have been calculated. The 0/C abundance and 
even more so Si/C increase to highly supersolar values 
if very massive stars leading to pair creation SN form. 
How far measurements of such abundance ratios, e.g. in 
the Lyman-a forest, can be used to constrain the IMF 
of the first stellar objects, remains currently unclear, 
as also discussed by Abia et al. (2001). 

The present model set should also be useful to ex- 
plore the feasibility of future direct observations of Pop III 
"galaxies" and to optimise search strategies for such ob- 
jects. We have started undertaking such studies (Schaerer 
& Pello 2001, Pello & Schaerer 2001). 

Tumlinson et al. (2001) speculate that some emission 
line objects found in deep Lyman-a surveys could actu- 
ally be He 11 A 1640 emitters such as the metal-free galaxies 
discussed here. Recently Oh et al. (2001a) have estimated 
the number of Pop III sources and mini quasars emitting 
strong He 11 lines. Although still speculative, their calcu- 
lations predict sufficient sources for successful detections 
with deep, large field observations foreseeable with the 
Next Generation Space Telescope or future instruments 
on ground-based 10m class telescopes. Pilot studies now 
also start to address the question of dust formation and 
obscuration in the early Universe (Todini & Ferrara 2001), 
which remain a concern for observations. In any case the 
upcoming decade should bring a great wealth of new in- 
formation on the early Universe, and possibly already the 
first in situ detection of the long sought Population III 
objects. 
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